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Recently we succeeded to make a reliable EFT prediction in a totally parameter-free 
manner for the S factors for the solar pp and hep processes, p + p — > d + e + + u e and 
3 He + p — > 4 He + e + + z/ e (I]. The strategy used in there is to embed a highly sophisticated 
standard nuclear physics approach (SNPA) exploiting realistic potentials |2l3j into an EFT 
framework, that we refer to as EFT*. Up to next-next-to-leading order in chiral expansion, 
it turnes out that there is effectively only one counter-term relevant to this process, the 
coefficient of which - (d R ) - has been renormalized to reproduce the experimental value 
of the tritium-beta decay. Our study has also led to very accurate EFT calculations on 
two-body weak processes that also receive contributions from the d R term, /i — d capture 
rate [3], and v — d scattering cross section [3]. 

1. Introduction 

In making theoretical predictions on electroweak processes of nuclear systems at low 
energy, we have basically two distinct approaches: One is the traditional model-dependent 
approach based on accurate nuclear potentials and the other is effective field theory (EFT) 
approach. The former, that we refer to as "standard nuclear physics approach (SNPA)", 
has so far scored an enormous success in wide areas of nuclear physics, achieving in some 
cases an accuracy that defies the existing experimental precision. This suggests that all of 
the essential ingredients of QCD - believed to be the mother theory of strong interactions - 
have been correctly encoded in SNPA. It however suffers from model-dependence and lack 
of systematic expansion scheme. In EFT approach, processes are described consistently 
and systematically, and free from the mentioned problems. Respecting these nice aspects, 
intensive works based on EFT are recently being made with impressive successes. For 
complicated systems with A = 3, 4, • • •, it is technically very hard to apply EFT directly 
and not much progress has been made up to date, although intensive efforts are being 
made in this direction. Now our question is, can we construct an EFT that exploits 
phenomenological but accurate SNPA wave functions, taking advantage of the merits of 
these two processes ? Our answer on this question is yes, and we have recently developed 
such a formalism which we call "more effective effective field theory (MEEFT)" or simply 
EFT* |Uj7|8|9| . In EFT*, it is possible to have the consistency and systematic aspect of 
EFT with the wave functions obtained accurately in SNPA. In this talk, I wish to report 
our developments on EFT*. 
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To be concrete, we shall discuss the following two solar nuclear fusion processes: 



Both figure importantly in the solar neutrino problems. Since the thermal energy of the 
interior of the Sun is of the order of keV, and since no experimental data is available 
for such low-energy regimes, one must rely on theory for determining the astrophysical 
S'-factors of the solar nuclear processes. Here we concentrate on the threshold ^-factor, 
S(0), for the reactions (JTJ and ©. The necessity of a very accurate estimate of the 
threshold S'-factor for the pp process, S pp (0), comes from the fact that pp fusion essentially 
governs the solar burning rate and the vast majority of the solar neutrinos come from 
this reaction. Meanwhile, the hep process is important in a different context. The hep 
reaction can produce the highest-energy solar neutrinos with their spectrum extending 
beyond the maximum energy of the 8 B neutrinos. Therefore, even though the flux of the 
hep neutrinos is small, there can be, at some level, a significant distortion of the higher 
end of the 8 B neutrino spectrum due to the hep neutrinos. This change can influence 
the interpretation of the results of a recent Super-Kamiokande experiment that have 
generated many controversies related to neutrino oscillations |lUlllj . To address these 
issues quantitatively, a reliable estimate of S/i ep (0) is indispensable. 

Before going further, I would like to explain the difficulty in making a reliable esti- 
mation of S'hep(O). First of all, the leading one-body contribution for the hep process is 
strongly suppressed due to the pseudo-orthogonality between initial and final wave func- 
tions. Secondly the main two-body (2B) corrections to the "leading" IB term tend to 
come with the opposite sign causing a large cancellation. The 2B terms therefore need to 
be calculated with great precision, which is a highly non-trivial task. Indeed, an accurate 
evaluation of the hep rate has been a long-standing challenge in nuclear physics ^U]- The 
degree of this difficulty may be appreciated by noting that theoretical estimates of the 
hep S'-factor have varied by orders of magnitude in the literature. 

2. Theory 

The primary amplitudes for both the pp and hep processes are of the Gamow- Teller 
(GT) type. Since the single-particle GT operator is well known at low energy, a major 
theoretical task is the accurate estimation of the meson-exchange current (MEC) contri- 
butions. In getting the current operators, we rely on the heavy-baryon chiral perturbation 
theory (HB^PT) with the Weinberg's power counting rule, which is a well-studied EFT 
that has been proven to be quite powerful and successful in describing low-energy nu- 
clear systems. In HB^PT, we have pions and nucleons as pertinent degrees of freedom, 
with all other massive degrees of freedom integrated out. The expansion parameter in 
HB^PT is Q/A x , where Q stands for the pion mass and/or the typical momentum scale 
of the system, and A x ~ 1 GeV. In our studies, we have calculated the MEC of the 
GT operator (space part of the axial-vector currents) up to N 3 LO, where N^LO stands 
for the order of (Q/A x ) u compared to the leading order (LO) one-body operator. Up to 
this order, MEC consists of one-pion-ranged and contact (zero-ranged) two-body currents, 
A2B = A l7T + A s . Three-body and many-body operators appear only in higher orders. 



pp : 
hep : 



p + p-^d + e + + v e , 

p + 3 He -»• 4 He + e + + u e . 



(1) 
(2) 
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A l7r can be determined unambiguously, thanks to available accurate ttN scattering data. 
A s however contains one unknown parameter, the overall strength of the contact-term 
contribution, which we denote by d R , 

A s cc d R J>, - Ttf-»{vi x tr^irti), (3) 

i<j 

where A is the cutoff (which we introduce in the procedure of Fourier transformation from 
momentum space to coordinate space), and 5^ (r) is the smeared delta-function with the 
radius ~ 1/A. d R represents all the heavy degrees of freedom integrated out, and chiral 
symmetry (or any other symmetry) does not tell us the value of it. We can however 
fix d R by studying other processes which are sensitive on it. In other words, we adjust 
the value of d R to reproduce the experimental data of those other processes, a procedure 
which is the renormalization condition of d R . This procedure is closely analogous to the 
EFT approach to effective nuclear potential backed by renormalization group equations as 
explained in [T2]. The power of the approach is that the constant d R appears in tritium 
/3-decay, /i-capture on a deuteron, and u-d scattering, as well as in pp and hep processes 
and hence is completely fixed. Among them, accurate experimental data is available for 
the tritium-beta decay rate, Pg, which we use to fix d R . 

In EFT*, we evaluate the transition matrix elements by sandwiching the obtained cur- 
rent operators between phenomenological but accurate SNPA wave functions. Recalling 
the fact that, the short- distance behaviors can be substantially different from model to 
model even among modern potentials, one may worry about model-dependence. To ad- 
dress this question, it is important to observe that, as was also recently discussed in Refs. 
[12], the dependence occurs only at short-distance, and thus the model-dependence is lo- 
cal. And we know that local (short-distance) contribution is well described (or controlled) 
in terms of the counter-terms, whose generic form is J2v Ci/Vtfjjf (r). The ci^-term is the 
leading counter-term. It means that, for a fixed A, the matrix element of (S^), may 
have substantial model-dependence (due to the model-dependent short-range behavior), 
which is, however, compensated by the model-dependence of the value of d R , to reproduce 
Fg. As a result, we can have the model-independent theory prediction for the total net 
amplitude. In other words, EFT* can be different from EFT only at short-distance, which 
is to be renormalized away order by order. 

3. Results 

To determine d R from P^, we calculate from the matrix elements of the current 
operators evaluated for accurate A=3 nuclear wave functions. We employ here the wave 
functions obtained in Refs. |3|13j using the correlated-hyperspherical-harmonics (CHH) 
method |14|15j . It is obviously important to maintain consistency between the treatments 
of the A=2, 3 and 4 systems. We shall use here the same Argonne v±s (AV18) potential jTHj 
for all these nuclei. For the A > 3 systems we add the Urbana-IX (AV18/UIX) three- 
nucleon potential ^7]. Furthermore, we apply the same regularization method to all the 
systems in order to control short-range physics in a consistent manner. 

In Tabled we have listed the value of d R and matrix elements of MEC, M.2B- The table 
indicates that, although the value of d R is sensitive to A, M.2B is amazingly stable against 
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Table 1 

The strength d R of the contact term and the two-body GT matrix element, M.2B, for the 
pp process calculated for representative values of A. 



A (MeV) 


d R 


M 2 b (fm) 


500 


1.00 ±0.07 


0.076 - 0.035 d R ~ 0.041 ± 0.002 


600 


1.78 ±0.08 


0.097 - 0.031 d R ~ 0.042 ± 0.002 


800 


3.90 ±0.10 


0.129 - 0.022 d R ~ 0.042 ± 0.002 



the variation of A within the stated range. The A-independence of the physical quantity 
Ai, which is in conformity with the tenet of EFT, is a crucial feature of the present EFT* 
result. Since the pp amplitude is dominated by the well-known one-body LO contribution, 
the resulting S factor at threshold, S pp (0), can be predicted with extreme accuracy, 

S pp (0) = 3.94 x (1 ± 0.0015 ± 0.0010) 10~ 25 MeV-b. (4) 

Here the first error is due to uncertainties in the input parameters in the one-body part, 
while the second error represents the uncertainties in the two-body part. 

In Table El we have listed results for the matrix element of the GT operators for the 
hep processes in arbitrary unit. We see from the table that the variation of the two-body 
GT amplitude (row labelled "2B-total") is ~10 % for the range of A under study. It is 
also noteworthy that the variation of the 2B amplitude as a function of A is reduced by a 
factor of ~7 by introducing the d R -term contributions; compare the third and fifth rows 
(labelled "2B (without rf R )" and "2B-total", respectively) in Table H As discussed, the A- 
dependence is amplified in total (IB + 2B) amplitude due to the substantial cancellation 
between IB and 2B contributions. The resulting 5-factor (adding the contributions from 
non-GT channels) reads 

S hep (0) = (8.6 ± 1.3) x 10- 20 keV-b , (5) 

where the "error" spans the range of the A-dependence for A=500-800 MeV. This result 
should be compared to that obtained by Marcucci et. al. Sh ep (0) = 9.64 x 10~ 20 keV-b. 



Table 2 

Values of the two-body GT amplitude (in arbitrary unit) for the hep process calculated 
as functions of the cutoff A. The individual contributions from the one-body (IB) and 
two-body (2B) operators are also listed. 



A (MeV) 


500 600 800 


IB 


-0.81 -0.81 -0.81 


2B (without d R ) 
2B (oc d R ) 


0.93 1.22 1.66 
-0.44 -0.70 -1.07 


2B-total 


0.49 0.52 0.59 
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4. Discussions 



The hen process, 3 He+n — > 4 He+7, shares many features with the hep process, including 
the suppression of the leading one-body contribution due to the pseudo-orthogonality of 
the wave functions. Especially, the hen process also contains two counter-terms, one in 
isoscalar and the one in isovector channel, which can be renormalized by the magnetic 
moments of 3 He and 3 H. We note that the MEC of the hen process starts from NLO 
coming from unsuppressed one-pion-exchange diagrams, while that of hep starts from 
N 3 LO. We are applying exactly the same strategy used in the hep process (work in progress 
with Y-H. Song). Accurate experimental data is available for the hen cross section, but so 
far no theoretical calculations have succeeded in explaining the data quantitatively. Thus 
applying the same EFT technique will provide a useful check of the validity of EFT* as 
a bona-fide EFT of QCD as discussed in Ref. |7|12| . 

This talk is based on the work done in collaboration with L.E. Marcucci, R. Schiavilla, 
M. Viviani, A. Kievsky, S. Rosati, K. Kubodera, D.-P. Min, M. Rho, S. Ando, F. Myhrer 
and Y.-H. Song, to all of whom I wish to express my sincere thanks. This work was 
supported by Korea Research Foundation Grant (KRF-2001-050-D00007). 
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